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1. INTRODUCTION 

Generator excitation control systems contain Automatic Voltage Regulators (AVR) for voltage 
regulation and conventional Power System Stabilizers (CPSS) for damping mechanical mode oscillations. 
The changes in operating conditions of PSS is challenge to update the controller parameters [1]. Therefore, 
the new studies seek to design advanced control techniques, which controllers adapt with the continuous 
change in operating points [2-4]. The conventional PID controller is common use in several of engineering 
applications. Due to the structure simplicity and easy parameter tuning, it is suitable for a certain operating 
point. In addition, its performance is good for linear and simple systems [5, 6]. Still, the behavior of PID 
control is linear and cannot deal with the high disturbance and high nonlinearity in complicated systems [5, 7, 
8]. The current research directed to use the Fractional Order PID (FOPID) control where it presents 
the nonlinear face of PID control [9-11]. In FOPID controller, two additional parameters (the fractional 
integral and derivative gains) will be supplementary to increase the flexibility and reliability of controller 
[12-14]. Therefore, the dynamic performance of FOPID controller is enhanced compared to 
the conventional PID controller [15-17]. 

At different operating points for a certain system, adaptation online was used self- tuning using for 
the system. In this case, the fuzzy logic calculations need a long time and addition efforts by try and error is 
performed to obtain normalizing gains selection [18]. So, this study resort to the MRAS to self-tuning 
the TSMFOPID online where it has simple structure, easy to implement and fast calculations [19, 20]. 
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The model-reference adaptive system (MRAS) presents one of the best adaptive control techniques. It may be 
regarded as an adaptive servo system in which the desired performance is expressed in terms of a reference 
model, which gives the desired response to a command signal [20]. It forces the overall system to follow 
the behavior of preselected model reference. The preselected model can be first or second order system 
according to the point of view the designer and complicated degree of the system [19]. 

In this Study, Takaji- Sugeno Fuzzy Type 3 (TS-Fuzzy Type 3) is designed by 21 rule triangular 
memberships. In the TSMFOPID design, the power of S operator for FOPID is lied between zero < À > 2 for 
integral fraction order and zero < u > 2 for derivative. Let the FOPID operated via the Ninteger Toolbox with 
internally unknown five parameters (kp, ki, kd, A and u) be named as (TBFOPID). While, a modified FOPID 
which has externally unknown five parameters (kp, ki, kd, A and p) constructed designed by TS fuzzy is 
named as TSMFOPID [21]. The performance of the TSMFOPID based on model reference adaptive system 
controller can be improved using different types of optimization technique Harmony Search (HS) [22, 23]. 
This paper presents, a new combination between modified FOPID controller based on TS technique 
(TSMFOPID) and Model reference as a tuner to design a new adaptively output feedback controller for PSS. 


2. POWER SYSTEM STABILIZER MODEL 

A single machine-infinite bus system whose linearized incremental model containing the voltage 
regulator and exciter can be demonstrated by the block diagram as shown in Figure 1. The parameters of 
the system are given in Table 1 [1, 24]. A number of cases are done that cover different operating points 
(normal, heavy and light load conditions) and parameters variation in the presence of a severe disturbance. 
These cases are applied to the system with optimal FOPID and self-tuning by Model Reference control. 
The constants K = [K1, ..., K6] in the normal, heavy and light loads are given in Table 2. The parameters of 
the system are given in Table 2 [2, 25]. 


Synchronous Machine 





Vo|tage Regulator 
and Exciter 





Figure. 1. Linearized incremental model of synchronous machine with an exciter and stabilizer. 


Table 1. The parameters of the system 
Parameter Value Parameter Value 





Ka 400 D 0 
Ta 0.05 Efdmax 7.3 
Tdo 5.9 Efdmin -7.3 
M 4.14 umax 0.12 
Kf 0.025 umin -0.12 
Tf 1 0 
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Table 2. Operating conditions for K_1 to K_6 
OP1=[1.0, 1.0] OP2=[1.3,0.9] OP3=[0.8, 1.2] 





K 





Normal load Heavy load Light load 
Kl 1.0753 0.6234 1.4076 
K2 1.2581 1.2813 1.1984 
K3 0.3071 0.3071 0.3071 
K4 1.7131 1.7123 1.6461 
K5 -0.0476 -0.2091 0.0742 
K6 0.4972 0.4565 0.5488 





3. CONTROL TECHNIQUES 

This section shows the proposed control techniques. The first technique is the Modified FOPID 
controller based on Harmony search. The second technique is the self-tuning of modified FOPID controller 
based on optimal MRAS. 


3.1. The FOPID Control 

The Toolbox FOPID (TBFOPID) is usually used to simulate the FOPID control. It has five 
parameters internally selected by designer in one closed block as shown in Figure.2. To make adapt 
the FOPID control online a Takaji- Sugeno Fuzzy (type-3 fuzzy) technique is developed FOPID has external 
five terminals parameters as shown in Figure 3 [21, 26]. 
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Figure 2. The block diagram with internal five Figure 3. The block diagram with the external five 
unknown parameters kp, ki, kd, A and u parameters to be suitable for model reference self- 
tuning 


The design steps of TSMFOI and TSMFOD can be summarized as follows [21]. The first step, let 
the input TS fuzzy membership functions for the fractional orders of the integral and derivative (A and u) 
values are select to be 21 triangular member ships functions. The universe of discourse values are equally 
distributed over the range [0, 2] and have their middle vertices placed at the points (0, 0.1, 0.2, 0.3, 0.4, 0.5, 
0.6, 0.7, 0.8, 0.9, 1, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2} .The membership selected by 21 triangular 
member ships as shown in Figure 4. The block of TSMFOD or TSMFOI and TSMFOPID 21 rules represents 
in Figure. 5. 
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Figure 4. Input membership of the variables of à and u 
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Figure 5. TSMFOPID controller with 21 memberships 


The second step, recognize the TS-Fuzzy formula for the fractional orders of integral and derivative 
parameters (A and u) as shown in Figure. 6. If the input is à the block diagram represents TSMFOI while if 
the input is u the block diagram represents TSMFOD. 


MM esght Calculabons 





Fractional zx 


Figure 6. Block diagram the fractional orders of the integral and derivative. 


The final outputs of the fuzzy systems that inferred for the TSMFOD or TSMFOI Implemented 
using Centroid for the defuzzification method [21]: 


Da Wieth 


J O tp T 
da, Wa; 


out, — 
Zu; Wy; 


; (1) 


Where: Ai wie (0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1, 1.1, 1.2, 1.3,1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 
2 for 21 rules 
Wh is the weight of ài, 
Wui is the weight of ui, 
F% is the output of TBFOPI whose A value is 4 
Fui is the output of TBFOPD whose u value is ui [12]. 
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3.2. Harmony search optimization technique 

The challenge point in the PID and FOPID controllers are selecting the appropriate parameters for a 
certain controlled plant. There are several methods to find the parameters of FOPID controller for example, 
try and error and Ziegler-Nichols method but, most of these techniques are rough roads. In this paper, the 
harmony search optimization technique will be used to obtain the optimal values of FOPID controller 
parameters according to the objective function as shown in (2) [22]. 


1 
— (1-e78)(My+ess)+e78 (ts—ty) 


f (2) 


Where e,, is the steady state error, M, is the overshoot of system response, t, is the settling time 
and t, 1s the rise time. Also, this objective function is able to compromise the designer requirements using 
the weighting parameter value (B). The parameter is set larger than 0.7 to reduce over shoot and steady state 
error. If this parameter is adjusting smaller than 0.7 the rise time and settling time will be reduced. Harmony 
search (HS) was suggested by Zong Woo Geem in 2001 [27]. It is well known that HS is a phenomenon- 
mimicking algorithm inspired by the improvisation process of musicians [28]. The initial population of 
Harmony Memory (HM) is chosen randomly. HM consists of Harmony Memory Solution (HMS) vectors. 
Table 3 shows the obtained parameters of TSMFOPID controller based on harmony search optimization 
technique. 


Table 3. TSMFOPID parameters. 
TSMFOPID parameters Kp Kd vd ki vi 
Parameters values 9.5603 5.3506 0.23714 2.5926 0.922 


Both of conventional toolbox of FOPID and the Takaji-Sugeno (TS) modified FOPID (TSMFOPID) 
have the same response through the simulation results at different operating conditions. In addition, it is 
provided in [21]. 


Kyu). Kiaz Kaa3z vday vias 
Key Kiez Ka(23; Vd(24); vigs) 


Kp (ums, 1) Ki (ums,2)Ka (us) vd (HMS,4) Vl (HMS,5) 


3.3. The self-tuning TSMFOPID based on model reference technique 

In this paper, the modified FOPID control parameters will be adjusted on-line using the model 
reference technique. The Model Reference Adaptive Control (MRAC) considers high-effectiveness adaptive 
controller [19]. It works as an adaptive servo system in which the wanted performance is described in form of 
a reference model. Figure. 7 demonstrates the main construction of self-tuning modified FOPID based on 
model reference technique. The details of Model Reference Adaptive System derivation are given in [29]. 






Adjustment 
Mechanism 


awe 


TSMFOPID Power System 
Uc Controller Stabilizer 





Figure 7. The overall system with self-tuning TSMFOPID based on model reference technique 


Model reference self-tuning fractional order PID control based on ... (M. A. Abdel Ghany) 


1338 O ISSN: 2088-8694 


The transfer function of FOPID control can be described as follows [29]. 
u(s) 
25 —k,yc k= A + kgs" (4) 
e=u.-y (5) 


Assume that the plant can be simplified to a first order system as obvious in the following equation. 





AS PME: 
u(s)  ‘Tst1 (6) 


Where k and T are unknown parameters. Also, assume that the model reference takes a form first 
order system as the following relationship. 


Ym(S) ES km (7) 


Uc(S) TmSt+1 


Where kn and Tm are selected by designer. 
From eduatons [4- 6l can conclude that 





= Ts+1 n t kis så + kas“) (uc — y) (8) 
yields kkptkkigtkkas" kkptkkigtkkas" 
———À — E Ey, e — — 
E E Ts+1 
kk, + kki. + tkkgs!\ kk, + kki. 2 + kkas” 
juge oe ee ee 
TS da Is i 
Ts 41-4 kk, + kk;. at kkas" kk, + kki. pt kkas" 
AME MM UNI po ou 
IS a Ts+1 
kkp+tkki- Z +tkkqas” 
y= ——— 12€ (9) 
ST LEkKy E kki-- j'kkas" 
= y — Ym (10) 
kky-kk;—-kkgs i 
= Ó——ÓÀ c (11) 
Ts+1+kkp+kki-j+kkąs” — Tms+1 
Bem _ " E K(kkptkkiGtkkgs") js 
ðkp Ts kkyFkkas!-kki +1 (rss aste n). c 
The (12) can be rewritten 
Beri K(Ts+kkp+kkj Gtkkgs!+1—kkp—kkiG—kkgs") : T 
ðkp (Ts+kkptkkit+kkast+1) 
ðe k(Ts+1) 
ET = E Uc (14) 
p (Ts+kkp+kki-ztkkast+1) 
dem k(Ts+1) 
ALL N 1 
ðkp rrr runs]? do 
From (4) and (6) 
dem k?e 
o —«— € — 1 
ðkp irn) ule) 
To achieve the desired performance, the following condition must be hold. 
Ts + kk, + kky. a+ kkas” +1 = Tus +1 (17) 
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dem _  k*e 
Ok, ^ TuS (9) 
From the MIT rule can obtain the following relationship 
ii M HE (19) 
d rn Cm TS 
dkp _ €m.€ 
a lA ml (21) 
yi — y. k* a (21) 
kp)new = ers dt + k,(0) (22) 
Where k, (0) is the initial value of proportional gain kp, By the same steps. 
dki 
ki)new = dt dt + kj(0) (23) 
Where k;(0) is the initial value of proportional gain k;. 
dkq 
ka)new = J =p at + ka(0) (24) 
Where k,,(0) is the initial value of derivative gain k4. 
dem _ kkiin(s)| Cete phas) E (25) 
OA så (Tstkkptkkyy+kkgst+1) (TstkkptkkiGtkkgst+1) : 
kKy-kkj--kkgqs" —l5—kkg—kkj-;-kkgqs!—1 
m = = (s) | p sÀ d - D sÀ d u (26) 
(Tstkkptkkjqtkkas"+1) 
Ton = = (s) | aos ; (27) 
(Ts+kkptkki-jtkkast+1) 
dem m kk;jln (s) —(Ts+1) 
0A sî dra y ee 
Also, from (4) and (6) 
dem _ — k^kjn(s) e 
dA — s^ rp] Ce) 
dem _ _ k^kjn(s) e 
QA — s^  '"Tgs*1 (30) 
dA — k? kjln (s) e 
at DUm sà 'Tas41 (31) 
dA €q.€ 
dec C Tms+1 (32) 
k?kj(0)1 
ya = y. T. = y. k; (0). In (s) (33) 
da 
A)new = p dt + A(0) (34) 
1 
TA kka. st.1n (5) kka. s".]n (s) (kk, um kkas”) 
ðu — 1 7 1 ze 
Ts + kkp + kk;-y b kkas" +1 — (Ts + kk, + kki -z + kkas" + 1) 

36 kkqstIn (s)(Ts-kky kk; 5 kkaqs!--1-kky kk; 5 kkas") 
— = | --A St? |, (35) 
x (Tstkkp+kkjGtkkgs+1) 
dem = kkq.s" 1n (s)(Ts+1) zu, (36) 
ou (TstkkptkkjGtkkgst+1) 
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Oem _ 
T 


kkq.s" Jn (s)(Ts+1) 
Ou — lcs tkk;l kk si-e1)(kk tkk; kk st) y (21) 
p UA d D LA d 


Also, from (4) and (6) 





dem _ k*kq.s" In (s).e 

ðu re (20) 

ðem _ k*kg.s"In (s).e 

Ou 75 a 2) 

du _ k?kg.s" In (s).e 

ae Vem Tmst1 (40) 

du _— €g.€ 

at l5 Imst (41) 

ys = y. k?. kq(0). S), In(s) = y3.kq(0). In (s) (42) 
d 

I)new = JG, dt + u(0) (43) 


4. SIMULATION RESULTS 

This section demonstrates the simulation results of fixed structure TSMFOPID and self-tuning 
TSMFOPID based on MRAS applied to PSS with different operating points (heavy and light parameters) 
through several types of disturbances. 


Case 1: Fixed structure TSMFOPID performance at different operating condition. 

The mechanical torque Tm and Vref increases suddenly with step value 5% in case heavy, light and 
normal parameters values. The results are demonstrated in Figure. 8. It is clear that fixed structure 
TSMFOPID response cannot adapt the changes in operating conditions. So, the self-tuning becomes essential 
to obtain high performance through several operating conditions and disturbances. 
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Figure 8. The system dynamic response with step change 0.05 


Case 2: Comparison between the fixed parameters TSMFOPID and self-tuning for TSMFOPID at light load 
condition. 

In case of light parameters values, the mechanical torque Tm and Vref exposes to step change 5% 
high. The results are demonstrated in Figure. 9. It is noted that self-tuning TSMFOPID (MRAS) response has 
low fluctuations, small overshoot and it reach to reference point in small time. 
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Figure 9. The system dynamic responses at light condition 


Case 3: Comparison between the fixed structure TSMFOPID and self-tuning TSMFOPID (MRAS) 
performance at heavy condition. 

The proposed controllers were investigated by comparing the dynamic responses of the PSS at step 
disturbance 5% in the mechanical torque ATm and AVref. Figure. 10 illustrates the results of this case. It is 
obvious that the dynamic response of self-tuning TSMFOPID (MRAS) has a good performance compared to 
fixed parameters TSMFOPID controller where it has smaller overshoot and small settling time. 
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Figure 10. System dynamic responses with reference tracking. 


Case 4: The performance of the fixed structure TSMFOPID and self-tuning TSMFOPID during the 
parameters variations and heavy load condition. 

To investigate the robustness of proposed controllers, the inertia coefficient increased to become 
M=1.5 of normal value and the disturbance for ATm and AVref represented by a 0.05 step change from zero 
to 2 seconds, then, decreased by 0.03 from 2 seconds to 4 seconds and finally decreased by 0.01 as shown in 
Figure. 11. Figure. 12 shows the system responses driven by self-tuning TSMFOPID (MRAS) and fixed 
parameters TSMFOPID controller when. It is clearly seen that the self-tuning TSMFOPID (MRAS) 
overcomes these variations and give good response with a small settling time, thus indicating 
the effectiveness of the self-tuning TSMFOPID (MRAS) over a wide range of parameter variation and 
change of operating conditions 


Model reference self-tuning fractional order PID control based on ... (M. A. Abdel Ghany) 


1342 O ISSN: 2088-8694 





0.05 
0.045 | 
0.04 


0.035 + 


Disturbance 
= 
© 
t 


0.025 
0.02 


0.015 ; 











0.01 ' i ; 
0 l 2 3 4 5 6 
Time ( Seconds ) 


Figure 11. Variable disturbance for ATm and AVref 
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Figure 12. System dynamic responses with reference tracking 


5. CONCLUSION 

A novel scheme of self-tuning Takaji-Sugeno Modified Fractional Order PID (TSMFOPID) baas on 
the Model Reference adaptive system (MRAS) applied on Power System Stabilizer (PSS). TS Fuzzy 
technique is used to construct a (TSMFOPID). The objective of Model Reference Adaptive System (MRAS) 
tunes the five parameters of Takaji-Sugeno Modified FOPID controller online. Different operating points for 
PSS were implemented to investigate the robustness of proposed controllers. The harmony optimization 
technique used to obtain the optimal parameters of proposed controllers. The simulation results provide that 
Self-Tuning TSMFOPID based on (MRAS) have better performance than the fixed parameters TSMOFOPID 
Controller. 
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